Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1055 1619 (N40E07 L= 74)   332     38        1 (1029/1600-1087/1641)   0.0 E09-E04 N37-N41
      2  723 1537 (N29E29 L= 52)   748    240        2 ( 665/1461- 805/1615) -29.4 E36-E20 N24-N33
      3  417 1431 (N17E49 L= 32)  2511    499        6 ( 364/1258- 471/1571) -43.5 E59-E40 N07-N25
      4 1095 1508 (N32E02 L= 79)   367    120        1 (1081/1451-1108/1558) -74.1 E03-W00 N28-N35
      5  909 1372 (N20E12 L= 69)   205     80        0 ( 865/1358- 948/1384)  21.8 E14-E09 N19-N21
      6 1719 1251 (N23W47 L=127)  1183    250        1 (1646/1135-1788/1359) -47.3 W40-W53 N16-N29
      7 1377 1240 (N18W20 L=101)   271     92        0 (1361/1198-1393/1288)  90.0 W18-W20 N15-N21
      8  924 1224 (N11E09 L= 72)   197     75        0 ( 898/1194- 957/1258) -45.0 E10-E06 N09-N13
      9 1100 1225 (N14W02 L= 83)   204     79        0 (1057/1213-1145/1236)  14.0 W00-W04 N12-N14
     10 1646 1110 (N13W40 L=121)   490    118        1 (1632/1057-1662/1175) -81.9 W39-W41 N10-N17
     11  701 1158 (N04E22 L= 59)   255    121        0 ( 650/1146- 766/1174)  23.2 E24-E17 N03-N06
     12 1733 1101 (N14W48 L=129)   433     99        0 (1708/1065-1752/1139)  68.2 W45-W49 N11-N16
     13 1701 1018 (N09W46 L=127)   313    133        2 (1666/ 999-1754/1046)   0.0 W42-W50 N07-N10
     14  204 1006 (S11E60 L= 21)  1165    502        0 ( 104/ 979- 307/1028)   8.4 E76-E49 S13-S08
     15  772  942 (S08E14 L= 67)   224     98        0 ( 748/ 890- 811/ 986)  63.4 E15-E12 S10-S04
     16 1710  650 (S14W55 L=136)   546    111        1 (1669/ 623-1740/ 668)   8.1 W50-W57 S15-S12
     17 1148  425 (S36W19 L=100)   740    157        3 (1064/ 401-1221/ 448)  15.3 W12-W24 S38-S33
     18  471  410 (S50E41 L= 40)   509    122        0 ( 442/ 384- 509/ 433)  31.0 E45-E35 S52-S46
  

    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