Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  862 1608 (N32E23 L=241)   247     47        1 ( 837/1587- 888/1630) -26.6 E25-E21 N31-N33
      2 1577 1531 (N42W34 L=298)   726    150        0 (1504/1488-1643/1562)  -8.1 W26-W40 N39-N43
      3  318 1499 (N16E63 L=201)   471     69        0 ( 288/1469- 347/1526)  45.0 E65-E60 N13-N17
      4  913 1345 (N16E12 L=252)   703    234        3 ( 885/1255- 946/1452) -63.4 E16-E10 N10-N22
      5  367 1419 (N12E53 L=211)   326     56        1 ( 354/1385- 387/1448) -63.4 E54-E51 N09-N13
      6 1248 1383 (N25W08 L=272)   713    177        0 (1178/1316-1296/1433) -29.1 W02-W11 N21-N27
      7  483 1318 (N07E40 L=224)   763    183        1 ( 468/1235- 503/1414) -61.4 E43-E36 N02-N13
      8  827 1314 (N13E17 L=247)   279     88        0 ( 788/1286- 871/1329)   0.0 E19-E13 N11-N13
      9 1860 1063 (N16W61 L=325)   425     71        0 (1839/1029-1876/1102)  78.7 W58-W63 N13-N18
     10  696  931 (S12E18 L=246)  1181    350        2 ( 637/ 791- 786/1087)  78.7 E20-E14 S21-S01
     11  848  864 (S13E08 L=256)   207     38        0 ( 823/ 846- 872/ 881)   0.0 E09-E06 S14-S12
     12  300  795 (S27E49 L=215)   448     82        2 ( 255/ 774- 339/ 817)  16.7 E54-E44 S28-S25
     13 1166  757 (S14W13 L=277)   387    157        0 (1095/ 705-1224/ 795) -18.4 W08-W17 S16-S12
     14  924  711 (S21E01 L=264)   338    104        1 ( 876/ 694- 986/ 729)  26.6 E03-W03 S22-S19
     15 1159  610 (S23W16 L=281)   229     63        0 (1134/ 589-1194/ 629)  45.0 W14-W18 S24-S21
     16 1468  599 (S18W37 L=302)   293     45        0 (1450/ 566-1486/ 629) -63.4 W35-W39 S20-S16
     17  976  530 (S32W06 L=271)   319    100        0 ( 923/ 520-1028/ 541)  15.9 W02-W09 S33-S30
  

    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