Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1022 1653 (N41E12 L=123)   505    124        2 (1003/1609-1047/1689) -71.6 E14-E09 N37-N43
      2 1427 1560 (N41W22 L=157)   601    141        0 (1378/1519-1471/1600)  41.2 W17-W26 N37-N45
      3  405 1449 (N17E51 L= 84)   630    131        1 ( 356/1415- 445/1497)  63.4 E55-E48 N14-N20
      4 1330 1462 (N32W14 L=149)   328     50        1 (1310/1431-1352/1489) -63.4 W12-W15 N30-N34
      5  769 1344 (N16E21 L=114)   223    103        0 ( 712/1328- 806/1361)  26.6 E24-E18 N13-N17
      6  979 1249 (N13E06 L=129)   386    105        1 ( 922/1226-1019/1284)  38.7 E09-E03 N10-N15
      7  997 1148 (N07E03 L=132)   788    304        2 ( 881/1040-1092/1231) -32.7 E11-W03 N01-N10
      8  637  944 (S10E23 L=112)   333    137        0 ( 569/ 924- 689/ 974)  14.0 E27-E19 S11-S07
      9  994  854 (S11W01 L=136)   721    186        3 ( 966/ 767-1021/ 923) -71.6 E01-W03 S15-S06
     10  443  565 (S38E38 L= 97)   468    123        0 ( 402/ 539- 485/ 598)  32.0 E41-E33 S40-S35
  

    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