Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1281 1683 (N37W33 L=239)   298     38        0 (1262/1659-1308/1699)   0.0 W30-W35 N35-N38
      2  602 1595 (N44E22 L=183)   534    131        2 ( 550/1558- 668/1627)  11.3 E27-E16 N41-N45
      3 1678 1488 (N17W55 L=261)   355     85        0 (1638/1457-1710/1512) -45.0 W52-W57 N15-N19
      4  841 1447 (N29E05 L=201)  1554    372        5 ( 683/1372-1026/1505)   2.4 E16-W08 N26-N30
      5  445 1359 (N30E35 L=171)  1630    420        2 ( 285/1317- 649/1407) -16.9 E52-E19 N23-N35
      6 1907 1102 (S09W67 L=273)   910    110        2 (1878/1062-1935/1156)  56.3 W62-W72 S11-S05
      7 1848 1031 (S12W58 L=264)   498     72        1 (1811/ 998-1880/1062)   7.1 W54-W62 S13-S10
      8 1839  956 (S16W57 L=262)   451     44        2 (1816/ 929-1858/ 981) -45.0 W54-W59 S17-S14
      9 1761  807 (S24W49 L=254)  1278    272        1 (1681/ 720-1819/ 929)  37.6 W41-W54 S28-S17
     10 1496  880 (S15W27 L=233)   277     75        0 (1457/ 859-1531/ 899)  11.3 W24-W29 S16-S14
     11  973  854 (S08E06 L=200)   422    125        0 ( 935/ 803-1027/ 886)  26.6 E08-E02 S11-S06
     12 1649  802 (S22W38 L=244)   319     51        2 (1634/ 770-1667/ 842) -78.7 W36-W39 S24-S19
     13  538  545 (S20E41 L=165)  1458    461        2 ( 407/ 424- 659/ 679)  24.0 E55-E28 S26-S13
     14  271  621 (S12E62 L=144)   526     38        2 ( 253/ 596- 291/ 643)  18.4 E65-E59 S13-S10
     15  328  461 (S22E67 L=139)  3139    211        3 ( 277/ 364- 379/ 537)  42.7 E76-E59 S29-S17
     16  865  469 (S31E21 L=185)   253     36        0 ( 848/ 442- 881/ 497)  63.4 E22-E19 S32-S28
  

    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