Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1518 1685 (N28W56 L= 80)  2682    323        7 (1382/1638-1677/1737) -22.6 W43-W69 N21-N31
      2 1600 1524 (N17W51 L= 75)   622    155        0 (1531/1479-1678/1561) -48.8 W47-W55 N12-N20
      3  228 1221 (N32E54 L=329)  1559    274        4 ( 140/1111- 311/1291) -28.3 E71-E45 N23-N37
      4 1863 1168 (S09W61 L= 84)  1132    193        3 (1826/1101-1898/1228)  54.5 W55-W64 S12-S05
      5 1798 1158 (S08W53 L= 77)   961    132        3 (1766/1098-1833/1225)  66.8 W50-W57 S11-S04
      6 1691 1167 (S05W44 L= 68)   435    156        2 (1657/1093-1736/1210)  24.0 W39-W48 S07-S02
      7 1397 1123 (N00W23 L= 47)   434    152        1 (1342/1081-1461/1149)  -7.1 W19-W27 S02-N01
      8  248 1047 (N21E52 L=331)   347     86        1 ( 214/1021- 286/1073) -23.2 E55-E48 N18-N22
      9 1111  921 (S04W02 L= 26)  1886    481        3 ( 973/ 742-1233/1062) -73.3 E02-W06 S16-N04
     10 1420  986 (S08W22 L= 45)   331     78        1 (1379/ 952-1450/1018) -78.7 W19-W23 S10-S05
     11  524  965 (N11E31 L=353)   648    193        0 ( 439/ 942- 614/1007) -28.6 E36-E25 N07-N14
     12 1587  733 (S27W29 L= 52)  1396    259        3 (1496/ 691-1698/ 790) -23.6 W22-W38 S30-S21
     13 1381  626 (S28W12 L= 36)  1503    301        3 (1316/ 490-1420/ 764)  87.0 W08-W14 S37-S18
     14  703  649 (S11E28 L=356)   455    110        2 ( 656/ 615- 761/ 686) -50.2 E29-E24 S13-S07
     15  743  500 (S20E30 L=353)  1395    293        3 ( 655/ 422- 804/ 636)  40.2 E36-E23 S24-S13
     16  625  598 (S12E34 L=349)   491     98        1 ( 580/ 565- 674/ 628)   0.0 E37-E30 S13-S10
     17  536  404 (S22E50 L=334)  1692    256        3 ( 443/ 333- 625/ 471)   8.7 E65-E39 S25-S18
     18 1267  349 (S45E02 L= 21)  2570    424       10 (1153/ 261-1400/ 471)  20.6 E14-W10 S49-S38
  

    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