Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1810 1288 (N32W52 L=199)   458     49        0 (1783/1266-1834/1309) -11.3 W49-W54 N31-N33
      2  367 1268 (S01E48 L= 99)   309     95        1 ( 346/1227- 397/1301) -18.4 E51-E44 S02-N00
      3 1630 1185 (N21W36 L=183)   767    142        1 (1587/1118-1668/1248) -49.4 W31-W39 N16-N23
      4 1843 1199 (N27W55 L=203)   443     62        1 (1818/1180-1862/1218)  31.0 W52-W57 N25-N28
      5  672  890 (S17E19 L=129)   302    104        1 ( 639/ 857- 728/ 921)  -9.5 E20-E14 S18-S15
      6  408  853 (S25E37 L=110)  1078    360        0 ( 322/ 765- 532/ 904)  36.9 E44-E28 S32-S20
      7 1198  809 (S11W15 L=162)   276     82        0 (1151/ 791-1243/ 828)  31.0 W12-W17 S12-S08
      8  970  692 (S23W04 L=151)   320    104        0 ( 919/ 671-1014/ 724)  51.3 W01-W06 S24-S19
      9 1618  518 (S18W54 L=201)   492    151        3 (1590/ 490-1652/ 555)  -7.1 W50-W58 S20-S15
     10 1014  465 (S36W14 L=161)   580    152        1 ( 931/ 451-1088/ 482)  20.0 W07-W18 S38-S34
     11 1466  425 (S28W47 L=194)   988    151        0 (1382/ 392-1518/ 477)  53.1 W41-W49 S31-S23
  

    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