Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1416 1714 (N45W48 L=225)  1219    221        2 (1331/1666-1508/1764)   2.0 W38-W66 N42-N46
      2 1194 1464 (N31W17 L=194)   888    149        3 (1130/1419-1264/1517) -45.0 W13-W21 N27-N35
      3 1753 1435 (N20W63 L=240)   421     60        0 (1720/1420-1785/1449) -15.9 W59-W66 N18-N21
      4  358 1356 (N31E49 L=129)  1893    373        6 ( 215/1320- 529/1398)   5.0 E67-E32 N28-N33
      5 1775 1317 (N13W59 L=237)   547    128        0 (1719/1280-1835/1346) -28.6 W53-W65 N09-N15
      6 1651 1289 (N14W46 L=224)   347    113        0 (1618/1246-1696/1330) -56.3 W44-W49 N10-N16
      7  214 1198 (N21E62 L=116)  1409    191        0 ( 179/1123- 249/1288) -37.6 E70-E56 N16-N27
      8  563  957 (N05E30 L=148)   357    103        2 ( 534/ 909- 590/1012) -80.5 E31-E27 N02-N08
      9 1016  943 (N01E01 L=176)   241    107        0 (1003/ 895-1034/ 995)  63.4 E02-W00 S01-N03
     10 1698  846 (S14W44 L=221)   395    103        0 (1665/ 801-1722/ 894)  51.3 W41-W46 S15-S10
     11 1297  864 (S07W15 L=192)   244     76        1 (1262/ 833-1319/ 886)   0.0 W12-W16 S08-S05
     12  401  667 (S12E47 L=131)   296     34        2 ( 381/ 647- 424/ 690) -33.7 E48-E44 S12-S10
     13  925  546 (S23E12 L=166)   688    201        2 ( 876/ 470- 974/ 628)  63.4 E15-E07 S28-S18
     14  425  589 (S17E47 L=130)   687     90        1 ( 391/ 550- 456/ 623) -68.2 E49-E44 S19-S14
     15 1498  578 (S28W28 L=206)   278     95        0 (1483/ 530-1511/ 622)  90.0 W27-W29 S31-S25
     16 1119  424 (S34E00 L=177)  2741    535        6 ( 946/ 347-1257/ 592)  22.2 E12-W10 S40-S23
     17  518  434 (S28E46 L=132)   543     85        0 ( 480/ 408- 559/ 463) -38.7 E48-E42 S30-S25
     18 1365  272 (S53W22 L=199)  3293    467        7 (1189/ 208-1580/ 351)   3.8 W03-W49 S58-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