Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1166 1738 (N41E13 L=175)   579    131        0 (1131/1697-1197/1793) -33.7 E18-E09 N38-N44
      2  542 1678 (N20E51 L=137)   417     51        0 ( 514/1659- 572/1691)  26.6 E53-E48 N18-N21
      3  591 1571 (N14E42 L=146)   590    136        1 ( 573/1512- 615/1630) -58.0 E45-E38 N10-N18
      4 1325 1077 (N03W15 L=204)   479    157        2 (1265/1049-1406/1097)  18.4 W11-W20 N01-N06
      5  298 1043 (S22E48 L=140)   302     98        0 ( 258/1023- 335/1060)  26.6 E51-E45 S23-S20
      6 1478  832 (S05W31 L=220)   858    165        2 (1456/ 755-1501/ 918) -60.9 W27-W34 S09-N00
      7 1872  802 (N06W66 L=254)   514    161        1 (1836/ 788-1909/ 819)  15.3 W60-W71 N04-N07
      8 1210  761 (S17W18 L=206)  1021    236        4 (1093/ 726-1317/ 807)  42.3 W11-W22 S21-S11
      9  907  508 (S40W08 L=197)  1781    348        3 ( 742/ 475-1088/ 551)  15.6 E04-W21 S44-S35
     10 1489  444 (S25W49 L=237)  1410    291        2 (1357/ 418-1593/ 488)  29.4 W39-W55 S29-S19
  

    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