Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1162 1603 (N33W16 L=270)   612     82        1 (1121/1568-1203/1629) -36.9 W12-W18 N30-N35
      2 1415 1216 (N07W25 L=280)   590    123        2 (1352/1171-1466/1240) -15.9 W21-W28 N04-N08
      3 1809  842 (S18W55 L=309)  1213    195        3 (1769/ 758-1852/ 916)  77.5 W50-W58 S23-S13
      4  307  858 (S05E49 L=206)   875    108        2 ( 274/ 808- 347/ 908) -71.6 E52-E45 S07-S01
      5 1889  820 (S20W65 L=320)  1066    122        4 (1862/ 763-1910/ 873) -80.5 W61-W68 S23-S16
      6 1756  660 (S29W53 L=307)   951    206        2 (1658/ 619-1842/ 685)   5.7 W43-W64 S31-S27
      7  982  508 (S33E08 L=247)  1073    385        0 ( 811/ 447-1148/ 550)   0.0 E21-W03 S35-S30
      8 1381  450 (S41W22 L=276)  1306    364        1 (1234/ 392-1545/ 532)  15.6 W09-W34 S43-S36
      9  673  389 (S38E35 L=220)   775    177        3 ( 587/ 355- 757/ 416)   0.0 E43-E27 S39-S36
     10 1525  339 (S52W41 L=296)  1422    172        4 (1458/ 313-1608/ 368)  -2.9 W32-W54 S53-S48
  

    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