Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1262 1467 (N28W26 L=183)   501    122        3 (1222/1429-1313/1499)  14.0 W22-W30 N26-N29
      2  519 1477 (N45E29 L=128)   357     66        1 ( 495/1452- 558/1496)   9.5 E31-E25 N43-N46
      3  450 1283 (N33E35 L=122)  1372    377        0 ( 363/1155- 563/1398) -48.2 E43-E26 N23-N42
      4  196 1218 (N32E59 L= 98)   533    121        1 ( 173/1202- 229/1244)  -9.5 E61-E55 N31-N33
      5  239 1074 (N23E55 L=102)  1470    291        3 ( 187/ 947- 289/1163) -74.1 E60-E50 N14-N28
      6 1673 1104 (S03W43 L=200)   871    158        1 (1613/1075-1735/1136) -24.0 W38-W47 S06-N00
      7 1523 1048 (S02W31 L=188)   518    155        2 (1500/ 973-1564/1112)  50.2 W28-W34 S06-N00
      8  665 1070 (N16E21 L=136)   317    115        0 ( 629/1050- 726/1089)   0.0 E23-E17 N15-N17
      9  845  902 (N03E13 L=144)  1074    202        2 ( 778/ 826- 919/1002) -70.0 E15-E09 S01-N10
     10 1538  830 (S16W28 L=185)   464     62        0 (1524/ 787-1555/ 875)  80.5 W26-W29 S18-S12
     11  830  578 (S16E21 L=136)   720    158        3 ( 754/ 522- 881/ 612) -49.4 E25-E18 S20-S12
     12  682  277 (S35E45 L=112)   975    158        0 ( 646/ 236- 712/ 320)  28.6 E51-E40 S38-S31
  

    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