Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1044 1639 (N43E03 L=280)  3438    661        7 ( 739/1532-1298/1726) -12.4 E33-W17 N36-N47
      2 1565 1511 (N36W41 L=323)  3081    494        6 (1363/1434-1755/1550)   0.0 W22-W64 N30-N39
      3  921 1196 (N12E07 L=275)  1964    432        2 ( 889/ 992- 945/1390) -85.2 E10-E04 N00-N24
      4  279 1178 (N06E54 L=229)  2368    394        3 ( 232/1008- 307/1352) -73.3 E61-E50 S04-N16
      5 1633 1292 (N22W42 L=324)   392     89        0 (1613/1246-1653/1328) -71.6 W40-W42 N18-N24
      6 1695 1265 (N20W47 L=330)   588    137        3 (1655/1218-1720/1328)  45.0 W42-W49 N17-N24
      7  511 1285 (N15E36 L=247)   330    104        1 ( 464/1255- 563/1309)  15.9 E39-E32 N12-N16
      8  425 1206 (N09E41 L=241)   674    189        1 ( 402/1143- 468/1275) -76.0 E43-E37 N05-N13
      9  481  762 (S18E35 L=248)  1505    432        2 ( 433/ 568- 516/ 954)  76.5 E40-E32 S30-S05
     10  201  891 (S12E61 L=221)   786     80        1 ( 172/ 848- 230/ 921) -45.0 E65-E57 S14-S10
     11 1286  772 (S12W18 L=301)   233     19        1 (1268/ 749-1304/ 798)  90.0 W16-W19 S13-S10
     12 1049  734 (S16W03 L=286)   314    109        0 (1021/ 705-1076/ 775)  78.7 W01-W05 S17-S12
     13  918  623 (S24E05 L=278)   457    139        0 ( 853/ 595- 980/ 659)  26.6 E08-W00 S25-S20
     14 1757  577 (S23W63 L=345)  2583    203        6 (1667/ 528-1851/ 618)   4.4 W51-W78 S26-S19
     15 1554  468 (S31W46 L=329)  1365    166        2 (1488/ 434-1626/ 508)  33.7 W39-W51 S33-S27
     16 1334  397 (S37W29 L=312)   571    109        2 (1280/ 370-1390/ 430)  20.6 W25-W33 S39-S34
     17 1089  318 (S46W11 L=294)   415     63        1 (1053/ 300-1122/ 337)  31.0 W08-W14 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