Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  188 1320 (N38E61 L=261)  1023    141        1 ( 171/1280- 211/1363) -78.7 E63-E58 N35-N40
      2 1032 1276 (N22W06 L=329)  1339    324        4 ( 888/1227-1173/1323) -36.3 E02-W13 N16-N27
      3  565 1230 (N28E26 L=297)   596     89        3 ( 523/1186- 601/1266) -50.2 E29-E23 N24-N30
      4 1836 1197 (S02W62 L= 25)   768     95        3 (1816/1152-1861/1250)  78.7 W59-W64 S05-N00
      5  193 1066 (N22E60 L=263)  2652    423        3 ( 135/ 903- 279/1209) -65.8 E66-E53 N11-N31
      6 1257  941 (S03W12 L=335)  2472    726        4 (1106/ 757-1488/1114)  18.4 W01-W29 S12-N03
      7  722 1026 (N13E18 L=305)   283     83        0 ( 680/1009- 765/1051) -31.0 E20-E15 N11-N14
      8 1322  773 (S14W13 L=336)   247     68        0 (1304/ 737-1337/ 816)  63.4 W11-W14 S15-S11
      9  301  564 (S10E64 L=259)  3326    442        8 ( 223/ 506- 404/ 634) -27.9 E72-E54 S14-S04
     10  731  520 (S18E29 L=294)  1380    339        2 ( 675/ 407- 789/ 621)  40.6 E37-E23 S24-S12
     11  618  304 (S32E48 L=275)  1109    185        1 ( 553/ 281- 660/ 342)   4.8 E54-E42 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