Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  667 1683 (N58E20 L=  2)  1989    240        4 ( 584/1613- 776/1727)  -2.4 E30-E06 N52-N61
      2 1217 1639 (N41W30 L= 51)  1998    457        4 (1014/1605-1420/1667) -21.2 W14-W46 N34-N47
      3 1608 1516 (N24W54 L= 75)   329     53        0 (1588/1488-1633/1543) -90.0 W52-W55 N21-N25
      4  248 1345 (N37E62 L=319)   711     66        1 ( 216/1315- 283/1369)   5.7 E67-E57 N35-N38
      5  461 1310 (N33E36 L=345)   250     48        0 ( 437/1292- 487/1324)   0.0 E38-E33 N31-N33
      6 1151 1208 (N15W11 L= 32)   749    202        2 (1091/1133-1195/1277) -68.2 W08-W12 N10-N20
      7 1788 1027 (S09W53 L= 74)   340     88        0 (1778/ 988-1799/1069)  76.0 W51-W54 S10-S06
      8 1565  685 (S25W30 L= 52)  1382    186        3 (1524/ 606-1605/ 765) -63.4 W27-W33 S30-S18
  

    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