Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  831 1847 (N44E49 L=323)  1076    158        1 ( 776/1802- 874/1885)  -6.3 E59-E41 N40-N46
      2  422 1605 (N16E61 L=311)  1079    160        0 ( 361/1543- 460/1642) -10.8 E73-E52 N12-N18
      3  954 1578 (N27E20 L=352)   610    167        1 ( 922/1513- 993/1633) -77.5 E22-E18 N22-N31
      4  696  924 (S17E17 L=355)   375    125        1 ( 647/ 901- 742/ 968)  -8.1 E20-E13 S18-S15
      5  619  904 (S20E22 L=350)   280     84        1 ( 599/ 879- 650/ 944) -45.0 E23-E19 S22-S18
      6  826  694 (S28E03 L=  9)   707    279        2 ( 691/ 662- 954/ 725)  33.7 E12-W04 S32-S22
      7 1467  401 (S27W52 L= 64)   891    106        2 (1434/ 366-1509/ 445) -59.0 W48-W56 S29-S24
      8 1311  311 (S38W48 L= 60)  1793    212        6 (1238/ 271-1403/ 360)  63.4 W38-W54 S42-S32
      9 1173  300 (S43W39 L= 51)   436     52        0 (1147/ 272-1196/ 328) -20.6 W34-W42 S44-S41
  

    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