Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1305 1558 (N33W01 L= 27)   308     85        1 (1267/1538-1340/1581)   0.0 E02-W03 N31-N34
      2 1510 1245 (N20W22 L= 49)   546    184        0 (1449/1150-1554/1313) -28.6 W16-W27 N15-N22
      3  195 1145 (S19E60 L=326)   610    148        3 ( 179/1082- 214/1218) -77.5 E63-E57 S23-S14
      4 1849 1023 (N19W54 L= 81)  1505    301        7 (1825/ 925-1880/1147) -56.3 W50-W58 N13-N26
      5  277 1074 (S22E50 L=336)   374     92        1 ( 245/1043- 293/1117)  68.2 E52-E48 S24-S19
      6 1185 1003 (S03W09 L= 36)   183     66        1 (1167/ 963-1204/1034) -76.0 W08-W10 S05-S01
      7  678  757 (S31E13 L= 13)  1293    477        7 ( 481/ 703- 887/ 830)  11.7 E29-W00 S33-S25
      8 1482  668 (S12W38 L= 64)   280     88        1 (1454/ 619-1507/ 712) -81.9 W36-W38 S15-S08
      9 1362  525 (S23W37 L= 63)   460     59        2 (1321/ 496-1402/ 551)  59.0 W34-W39 S26-S21
  

    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