Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1049 1729 (N51W26 L= 59)   305     66        0 (1026/1704-1080/1748)   7.1 W22-W30 N49-N51
      2  474 1401 (N41E33 L=360)   309     66        1 ( 444/1379- 500/1422)  21.8 E36-E30 N40-N42
      3  193 1093 (N25E63 L=330)   966    157        4 ( 157/1014- 220/1154) -45.0 E69-E59 N20-N29
      4  597  964 (N12E27 L=  6)  1420    379        6 ( 541/ 843- 695/1136) -75.3 E30-E22 N03-N22
      5  243  926 (N15E58 L=335)   597    127        2 ( 219/ 871- 266/ 983)  40.6 E61-E54 N11-N17
      6  299  860 (N10E53 L=340)   500    118        0 ( 256/ 821- 340/ 894)  16.7 E58-E48 N08-N11
      7 1280  562 (S27W06 L= 38)   880    250        2 (1202/ 499-1327/ 680)  42.0 W00-W10 S29-S19
      8  864  484 (S22E23 L= 10)   530    149        1 ( 828/ 418- 905/ 540)  37.9 E27-E18 S26-S19
      9  588  291 (S31E56 L=337)   543     65        1 ( 564/ 263- 613/ 316)  21.8 E60-E50 S33-S29
  

    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