Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  543 1544 (N20E44 L= 42)   591    151        1 ( 499/1522- 607/1566)   6.3 E48-E39 N18-N22
      2 1348 1174 (N14W17 L=103)   377    153        1 (1293/1118-1408/1239) -29.1 W12-W21 N11-N16
      3  731  941 (S12E16 L= 70)   321    144        0 ( 710/ 870- 746/1012) -63.4 E18-E14 S16-S08
      4 1505  849 (S02W33 L=119)   426    124        1 (1477/ 795-1529/ 916) -81.9 W31-W34 S05-N02
      5 1878  878 (N07W67 L=153)   855    163        1 (1842/ 839-1903/ 908) -18.4 W61-W72 N05-N08
      6 1364  832 (S06W24 L=110)   215     56        1 (1351/ 802-1377/ 866) -56.3 W22-W25 S07-S04
      7  458  795 (S26E33 L= 53)   249    107        1 ( 422/ 783- 511/ 812)  15.9 E36-E29 S27-S24
      8  976  605 (S28W06 L= 92)   856    333        3 ( 801/ 550-1168/ 655)  36.9 E04-W16 S35-S20
      9 1630  537 (S18W54 L=140)   603    116        2 (1597/ 504-1665/ 576)  71.6 W51-W55 S20-S14
     10 1632  443 (S23W59 L=145)   607    181        2 (1580/ 421-1688/ 461)  21.8 W54-W64 S25-S21
     11 1461  254 (S40W58 L=144)   823    190        1 (1426/ 237-1504/ 273)  14.0 W54-W63 S41-S38
  

    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