Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1534 1606 (N23W55 L=342)  1248    368        3 (1408/1511-1658/1730) -84.3 W52-W58 N13-N33
      2  911 1556 (N40W10 L=297)   724    144        2 ( 850/1514- 969/1586) -45.0 W05-W13 N35-N42
      3 1043 1544 (N35W18 L=305)   483     83        3 (1007/1510-1076/1574)  11.3 W15-W20 N33-N37
      4 1226 1497 (N27W27 L=314)   683    272        0 (1106/1446-1343/1544) -52.4 W21-W32 N20-N33
      5  340 1519 (N56E50 L=237)  1151    119        2 ( 293/1485- 404/1540)   0.0 E61-E39 N54-N58
      6 1713 1401 (N06W57 L=344)   768    126        1 (1682/1350-1743/1459)  32.0 W51-W60 N04-N09
      7 1400 1355 (N14W32 L=319)   327     54        1 (1370/1332-1428/1379) -76.0 W29-W33 N11-N15
      8 1167 1228 (N14W14 L=301)   341    102        2 (1122/1198-1214/1259) -59.0 W11-W16 N11-N16
      9  203 1177 (N35E61 L=227)   817    111        2 ( 159/1138- 249/1217) -21.8 E69-E53 N31-N37
     10  771 1135 (N19E12 L=276)   281    100        0 ( 730/1094- 799/1169) -63.4 E13-E10 N16-N22
     11 1712 1074 (S11W44 L=331)   291     33        0 (1693/1052-1731/1097)  71.6 W42-W45 S12-S09
     12 1859  973 (S22W57 L=344)   429    105        0 (1821/ 954-1890/1000) -18.4 W52-W61 S23-S19
     13  757  923 (N08E18 L=270)   360    126        1 ( 709/ 900- 808/ 958)   9.5 E20-E13 N06-N08
     14 1612  722 (S28W29 L=316)   835    328        0 (1551/ 595-1702/ 841)  34.5 W21-W38 S34-S23
     15 1381  585 (S29W09 L=297)   971    213        3 (1306/ 521-1446/ 678)  30.3 W03-W15 S31-S24
     16 1714  584 (S42W38 L=326)  1248    104        2 (1682/ 535-1746/ 641)  76.0 W34-W42 S45-S37
     17  691  529 (S13E34 L=253)   273     69        1 ( 659/ 513- 728/ 545) -11.3 E36-E31 S15-S11
     18 1074  233 (S45E26 L=261)  1544    236        3 ( 964/ 183-1164/ 287)   5.0 E39-E16 S47-S43
  

    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