Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1110 1723 (N52W07 L=211)   390     70        1 (1076/1700-1152/1740)   7.1 W03-W11 N50-N53
      2  542 1519 (N34E40 L=164)  1794    396        2 ( 393/1448- 710/1605) -11.9 E61-E23 N29-N39
      3  723  688 (S18E20 L=184)   494    214        0 ( 625/ 666- 815/ 720)  -8.7 E26-E13 S19-S16
      4  384  586 (S27E49 L=154)  2331    323        4 ( 311/ 433- 457/ 692)  29.1 E77-E39 S39-S19
      5 1702  647 (S21W52 L=255)   452     41        3 (1679/ 625-1727/ 669)  71.6 W49-W55 S22-S19
      6  772  610 (S23E17 L=187)   441    162        0 ( 729/ 550- 838/ 648)  37.9 E20-E11 S27-S20
      7 1156  496 (S31W10 L=214)   281     84        1 (1110/ 479-1195/ 515)  31.0 W06-W12 S32-S29
      8 1491  384 (S40W42 L=246)  1851    259        2 (1397/ 336-1580/ 462)  45.0 W33-W50 S44-S34
      9 1206  309 (S47W18 L=221)  1164    205        3 (1126/ 285-1299/ 347)   3.6 W10-W26 S48-S43
  

    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