Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  834 1514 (N23E24 L=277)   283    122        0 ( 807/1457- 853/1560) -45.0 E26-E20 N20-N26
      2  536 1526 (N18E45 L=256)   283     80        0 ( 522/1493- 556/1554) -63.4 E46-E43 N15-N19
      3 1189 1452 (N27W01 L=302)   944    314        0 (1064/1360-1303/1508)  -9.5 E08-W10 N23-N30
      4 1371 1426 (N30W14 L=315)   304     98        1 (1337/1385-1406/1471) -33.7 W10-W16 N27-N32
      5 1869 1186 (N27W63 L=  4)  2247    225        2 (1827/1086-1898/1291) -90.0 W56-W68 N20-N34
      6  433 1067 (S11E39 L=263)   490     88        2 ( 382/1038- 466/1092)  33.7 E42-E36 S13-S09
      7 1754 1001 (N13W50 L=351)   681     98        0 (1733/ 948-1780/1042)  63.4 W47-W51 N09-N15
      8 1432  550 (S21W39 L=340)  2114    441        1 (1356/ 435-1548/ 657)  72.6 W34-W48 S28-S11
      9  460  557 (S44E34 L=267)   339     43        1 ( 447/ 529- 479/ 586)  78.7 E35-E31 S46-S41
     10  914  538 (S34W04 L=305)   644    212        1 ( 819/ 508-1004/ 566)  28.6 E02-W09 S37-S29
  

    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