Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1116 1680 (N34E12 L=307)   873    171        1 (1045/1625-1187/1714)  -4.8 E17-E05 N32-N36
      2  290 1395 (S00E58 L=261)   471     87        0 ( 263/1368- 329/1436) -12.5 E62-E53 S01-N01
      3 1378 1392 (N23W12 L=331)   231     81        0 (1366/1361-1394/1434) -53.1 W10-W13 N20-N24
      4 1233 1387 (N18W03 L=322)   396     85        0 (1186/1361-1283/1409)   9.5 W00-W06 N16-N19
      5 1670 1275 (N24W34 L=354)   408     85        0 (1646/1233-1701/1317) -38.7 W32-W37 N22-N26
      6 1876 1094 (N20W57 L= 16)  2425    304        3 (1841/ 983-1929/1252) -86.4 W52-W66 N14-N30
      7 1683 1097 (N14W38 L=357)   677    195        2 (1627/1031-1724/1152) -18.4 W33-W42 N10-N16
      8  609 1095 (S12E26 L=293)   202     68        1 ( 588/1065- 628/1126)  76.0 E26-E24 S14-S10
      9 1437 1051 (N04W23 L=342)   426    140        0 (1365/1026-1499/1087)  40.6 W18-W25 N01-N07
     10  538  874 (S26E27 L=292)   348     93        0 ( 498/ 854- 583/ 894)   9.5 E29-E23 S27-S25
     11 1826  747 (S00W61 L= 21)   694     87        1 (1786/ 730-1876/ 769)  10.3 W56-W68 S01-N01
     12  306  641 (S46E50 L=269)  2096    313        4 ( 176/ 593- 456/ 680)   0.0 E73-E32 S48-S43
     13 1733  563 (S12W62 L= 21)  1380    189        3 (1691/ 477-1775/ 645) -52.1 W57-W68 S17-S08
     14  592  546 (S46E19 L=301)   244     54        0 ( 573/ 515- 612/ 575) -36.9 E20-E16 S47-S44
     15 1658  490 (S19W60 L= 19)   531     85        1 (1641/ 438-1675/ 531) -24.0 W56-W65 S20-S16
     16 1486  468 (S25W47 L=  6)   560     98        0 (1439/ 450-1534/ 490)  14.0 W42-W51 S26-S23
     17 1473  315 (S34W59 L= 19)  2588    292        5 (1332/ 286-1610/ 359)  35.8 W48-W72 S39-S26
  

    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