Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1342 1672 (N32W43 L=102)   411    146        0 (1286/1638-1391/1702)   0.0 W37-W48 N31-N33
      2 1425 1494 (N20W39 L= 97)   320     82        1 (1405/1460-1445/1533)  18.4 W36-W42 N18-N21
      3  780 1197 (N22E09 L= 49)  1828    540        2 ( 554/1134-1035/1316) -41.3 E22-W04 N12-N35
      4 1762  975 (S18W45 L=103)  1901    454        2 (1715/ 779-1800/1197)  55.7 W37-W53 S28-S06
      5 1836 1119 (S13W56 L=114)   887    131        1 (1807/1069-1864/1159)  38.7 W52-W59 S15-S10
      6 1359 1128 (N03W21 L= 80)   333     58        3 (1332/1099-1382/1155) -76.0 W19-W22 N00-N04
      7 1542  717 (S26W23 L= 81)  1805    456        0 (1414/ 599-1601/ 955)  48.4 W12-W31 S32-S13
      8  870  773 (S04E15 L= 43)   286    102        0 ( 841/ 724- 918/ 825) -81.9 E16-E13 S07-N00
      9 1023  721 (S10E08 L= 50)   429    148        0 ( 948/ 689-1095/ 766)   0.0 E12-E03 S11-S09
     10 1093  632 (S18E07 L= 51)   275    128        0 (1067/ 572-1118/ 698)  45.0 E09-E03 S20-S14
     11  226  649 (N03E68 L=350)   524    188        0 ( 199/ 620- 253/ 685)  13.0 E75-E62 N01-N04
  

    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