Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  659 1775 (N58E23 L=154)   359     24        2 ( 638/1755- 681/1792)  -9.5 E25-E19 N55-N59
      2 1747 1574 (N24W70 L=246)   567     57        0 (1722/1557-1772/1592) -15.9 W65-W75 N22-N24
      3  237 1485 (N40E67 L=110)   599     54        1 ( 208/1463- 270/1500)  -9.5 E77-E59 N38-N41
      4  675 1452 (N31E18 L=159)   251     87        0 ( 649/1418- 707/1486) -51.3 E19-E15 N28-N33
      5 1260 1431 (N21W21 L=197)   221     69        0 (1240/1398-1283/1464) -78.7 W19-W21 N18-N23
      6 1604 1291 (N08W41 L=217)   326     80        1 (1580/1262-1636/1320)  18.4 W38-W44 N06-N09
      7 1460 1229 (N06W30 L=206)   218     76        0 (1443/1187-1480/1268)  45.0 W27-W31 N04-N08
      8 1795  899 (S18W53 L=229)   475    146        0 (1766/ 839-1824/ 969)  54.5 W49-W55 S20-S13
      9 1213  787 (S17W08 L=185)   941    417        1 (1155/ 630-1260/ 958)  78.1 W04-W10 S26-S06
     10 1664  881 (S17W40 L=217)   595     99        1 (1632/ 847-1714/ 912)   0.0 W37-W44 S18-S14
     11  571  614 (S18E36 L=141)   300     84        0 ( 552/ 574- 607/ 639)  31.0 E37-E32 S20-S16
     12 1592  407 (S49W43 L=219)   890    118        1 (1544/ 383-1644/ 441)  -7.6 W35-W50 S51-S45
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y