Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  613 1631 (N19E43 L=168)   463    153        1 ( 588/1572- 635/1684) -54.5 E47-E41 N15-N22
      2  484 1608 (N16E52 L=159)   842    173        2 ( 414/1580- 549/1641)  21.8 E58-E47 N13-N18
      3  788 1575 (N19E29 L=183)   465    172        0 ( 733/1514- 840/1634) -28.6 E34-E23 N16-N22
      4  684 1513 (N13E33 L=178)   624    231        1 ( 649/1421- 716/1589) -70.0 E38-E30 N07-N18
      5 1138 1518 (N24E06 L=205)   460    152        1 (1087/1451-1173/1577) -36.9 E10-E02 N20-N27
      6  404 1482 (N06E52 L=159)   357     74        0 ( 372/1464- 430/1511)  45.0 E54-E50 N04-N08
      7 1770 1442 (N40W46 L=257)   682     80        1 (1739/1406-1795/1478) -26.6 W41-W49 N37-N42
      8  739 1360 (N05E24 L=187)   255     69        0 ( 705/1341- 780/1380)   0.0 E26-E21 N04-N06
      9  176 1065 (S21E64 L=148)   609     82        1 ( 152/1025- 196/1101)  63.4 E67-E61 S23-S19
     10 1140 1068 (S02W05 L=217)   200     68        1 (1122/1036-1160/1097)  90.0 W04-W06 S04-N00
     11  181  861 (S34E64 L=148)   481    105        1 ( 162/ 810- 194/ 915)  45.0 E67-E60 S36-S30
  

    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