Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1088 1833 (N49E16 L=195)  1270    240        4 (1050/1768-1127/1898) -75.1 E19-E12 N43-N58
      2  368 1495 (N12E54 L=157)   738    158        1 ( 312/1437- 431/1534)  40.6 E58-E49 N07-N14
      3 1463 1405 (N25W21 L=232)   689    161        0 (1421/1347-1499/1455)  66.0 W18-W23 N19-N28
      4  370 1375 (N04E49 L=161)   602    101        0 ( 348/1316- 402/1430) -81.9 E51-E48 N00-N07
      5 1753 1289 (N24W44 L=255)  1216    159        1 (1708/1213-1791/1370) -60.9 W40-W48 N20-N29
      6  814  968 (S13E11 L=199)  1217    321        4 ( 699/ 833- 889/1113)  69.8 E16-E08 S22-S03
      7 1644 1038 (N06W38 L=249)   505    130        1 (1593/1001-1690/1089) -18.4 W33-W42 N04-N08
      8  335  927 (S23E45 L=166)   994    303        2 ( 240/ 834- 401/1015)  42.9 E53-E39 S29-S16
      9  436  770 (S31E36 L=175)  1399    419        1 ( 304/ 689- 574/ 891)  36.0 E48-E25 S37-S21
     10 1651  655 (S15W49 L=260)  4681   1070        7 (1261/ 588-1894/ 717)  23.3 W24-W78 S27-S05
     11  974  541 (S36W07 L=218)   254     80        0 ( 926/ 528-1015/ 553)  80.5 W03-W09 S37-S34
  

    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