Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  977 1709 (N51W21 L=290)   410     63        1 ( 954/1668- 995/1743)  35.5 W18-W25 N47-N53
      2 1379 1675 (N35W48 L=317)   398     94        1 (1331/1656-1421/1690) -33.7 W45-W51 N32-N37
      3 1662 1580 (N20W65 L=334)  1443    113        4 (1618/1533-1704/1626)  80.5 W61-W71 N17-N23
      4 1525 1546 (N23W50 L=319)   547    156        0 (1500/1462-1540/1610)  45.0 W45-W53 N18-N26
      5  452 1325 (N38E34 L=235)  1097    200        1 ( 364/1260- 529/1356) -29.4 E43-E27 N32-N41
      6  881 1275 (N25E02 L=267)   790    198        1 ( 780/1248- 969/1299) -14.0 E09-W03 N22-N27
      7 1070 1276 (N20W09 L=278)   254    104        0 (1033/1242-1117/1299) -59.0 W07-W11 N17-N22
      8 1452 1264 (N10W32 L=301)   213     38        0 (1431/1246-1476/1280)   0.0 W29-W33 N08-N10
      9 1285 1169 (N09W19 L=288)   651    151        1 (1224/1127-1336/1212)   6.3 W13-W22 N07-N11
     10 1059 1180 (N15W06 L=275)   194     32        0 (1037/1166-1085/1198) -33.7 W04-W07 N14-N16
     11 1682 1110 (S05W43 L=312)  1086    185        3 (1648/1025-1710/1190)  41.6 W38-W47 S09-N00
     12  741 1128 (N19E15 L=255)   224     69        0 ( 707/1108- 771/1151)   0.0 E16-E12 N18-N20
     13  432 1062 (N22E38 L=231)   931    169        2 ( 365/1002- 498/1113) -45.0 E43-E33 N17-N26
     14 1575  586 (S35W26 L=295)  5507    835        8 (1366/ 345-1789/ 928)  32.1 W06-W49 S48-S18
     15 1129  589 (S21E04 L=265)   515    217        1 (1071/ 523-1162/ 684)  45.0 E09-E01 S23-S15
     16  525  396 (S20E53 L=216)   452     73        0 ( 503/ 350- 541/ 419)  59.0 E55-E50 S23-S18
     17  978  341 (S34E23 L=246)   425     63        0 ( 946/ 314-1003/ 376)  23.2 E26-E19 S35-S31
     18 1283  296 (S48E04 L=265)   397     72        1 (1265/ 263-1309/ 327)  45.0 E06-W00 S51-S45
     19 1184  237 (S50E15 L=254)   638     90        3 (1140/ 214-1234/ 259) -20.6 E18-E09 S52-S47
  

    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