Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1272 1803 (N52W14 L=261)   669    101        3 (1248/1752-1293/1843)  84.3 W11-W17 N47-N57
      2  677 1684 (N35E34 L=214)   302     83        1 ( 645/1646- 702/1713) -29.7 E37-E30 N32-N37
      3  545 1659 (N32E44 L=204)   760    113        1 ( 518/1594- 565/1701) -63.4 E48-E40 N27-N35
      4 1472 1483 (N28W27 L=274)   444    128        1 (1427/1445-1523/1511)  15.9 W23-W30 N25-N30
      5 1645 1458 (N29W41 L=289)   292     44        2 (1628/1436-1671/1482) -90.0 W39-W43 N27-N30
      6 1721 1185 (N12W45 L=293)   423     59        0 (1688/1151-1744/1212) -26.6 W42-W47 N09-N13
      7 1678 1070 (N04W42 L=290)   288     63        1 (1658/1033-1700/1095)  76.0 W40-W43 N02-N06
      8 1708  642 (S20W54 L=302)   531     54        3 (1682/ 619-1734/ 669) -26.6 W50-W57 S21-S18
      9 1539  362 (S41W56 L=303)  1028     87        2 (1503/ 317-1566/ 415) -32.5 W50-W63 S43-S36
     10  568  310 (S57E43 L=204)  1140    152        4 ( 514/ 274- 625/ 351)  14.0 E56-E32 S60-S53
  

    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