Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1248 1398 (N20W24 L= 98)   436    225        2 (1223/1335-1271/1456)  35.5 W20-W27 N17-N23
      2 1814 1265 (S05W59 L=133)  2067    250        4 (1772/1168-1853/1385)  67.4 W55-W62 S10-N02
      3 1048 1192 (N15W06 L= 80)   431    171        0 ( 964/1173-1140/1225) -16.7 W01-W11 N12-N17
      4  647 1058 (N18E21 L= 53)   722    233        0 ( 550/1031- 767/1095) -23.2 E27-E13 N13-N21
      5 1583  842 (S20W28 L=102)  1525    434        4 (1518/ 665-1619/1074)  54.5 W20-W35 S28-S07
      6 1687 1019 (S13W40 L=114)   277     73        1 (1670/ 976-1703/1052)  53.1 W38-W41 S15-S11
      7  174  872 (N17E64 L= 10)   603    116        1 ( 160/ 825- 188/ 929)  56.3 E67-E61 N14-N20
      8  633  702 (S02E31 L= 43)   235     41        1 ( 621/ 671- 647/ 733)  56.3 E32-E29 S03-N00
      9 1333  568 (S29W05 L= 79)  1207    270        1 (1218/ 500-1403/ 678)   8.1 E02-W12 S32-S23
     10  762  600 (S11E27 L= 47)   275    109        1 ( 720/ 582- 817/ 628)   0.0 E29-E22 S11-S08
     11  876  481 (S21E25 L= 49)   710    263        0 ( 856/ 384- 899/ 597)  52.4 E28-E18 S27-S14
     12  540  423 (S16E49 L= 24)   524     97        1 ( 496/ 389- 586/ 448) -51.3 E52-E47 S19-S14
  

    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