Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  956 1830 (N64E08 L=  5)   621     55        0 ( 926/1808- 991/1850)  12.5 E12-E03 N61-N66
      2  958 1605 (N42E05 L=  9)   828    180        0 ( 895/1534-1016/1674) -42.3 E11-W00 N37-N48
      3 1177 1329 (N23W10 L= 24)   306    145        0 (1123/1284-1227/1372) -45.0 W07-W13 N19-N25
      4 1345 1318 (N22W22 L= 36)   248     85        0 (1298/1305-1378/1340)   9.5 W18-W24 N20-N23
      5  523 1270 (N19E34 L=339)   282     68        1 ( 491/1252- 568/1286)  -8.1 E37-E30 N17-N19
      6 1692  718 (S17W48 L= 62)  3475    566        5 (1589/ 556-1800/ 894)  69.1 W41-W57 S28-S06
      7  497  825 (S09E35 L=338)   237     26        1 ( 478/ 808- 518/ 844)  18.4 E36-E33 S09-S07
      8 1604  758 (S14W40 L= 53)   766    109        1 (1557/ 725-1652/ 804)  15.9 W36-W43 S16-S11
      9 1225  696 (S17W13 L= 26)   718    206        1 (1143/ 641-1317/ 744)  28.6 W07-W18 S20-S14
     10  640  493 (S31E29 L=344)  2528    346        6 ( 563/ 402- 715/ 588) -61.7 E34-E23 S38-S24
  

    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