Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1096 1688 (N34E14 L=307)   889    168        2 (1023/1637-1164/1721) -14.0 E19-E07 N32-N36
      2  352 1547 (N10E60 L=261)   372     82        3 ( 316/1515- 379/1579)  68.2 E62-E57 N07-N12
      3 1354 1403 (N22W10 L=331)   208     76        0 (1337/1369-1369/1442) -76.0 W08-W11 N20-N24
      4 1205 1399 (N18W01 L=322)   389    101        0 (1161/1373-1256/1423)   0.0 E01-W04 N16-N19
      5  318 1385 (S01E55 L=266)   676    108        0 ( 263/1369- 372/1401)   5.7 E60-E49 S01-N00
      6 1856 1235 (N28W54 L= 15)   535     66        0 (1835/1204-1880/1262)  45.0 W51-W57 N25-N31
      7 1669 1097 (N13W37 L=358)   601    153        1 (1646/1041-1708/1165) -39.8 W34-W40 N10-N17
      8 1877 1059 (N18W58 L= 18)  2366    254        7 (1830/ 995-1949/1165) -14.7 W50-W70 N14-N23
      9 1398 1059 (N04W20 L=341)   593    162        1 (1317/1030-1475/1096)  45.0 W15-W24 N00-N07
     10  516  884 (S26E29 L=292)   288     52        0 ( 485/ 863- 546/ 905) -11.3 E31-E26 S27-S24
     11 1814  756 (S00W60 L= 21)   617     75        0 (1775/ 742-1864/ 774)  10.3 W55-W66 S01-N01
     12  292  643 (S46E52 L=269)  2278    349        1 ( 165/ 597- 446/ 691)   0.0 E76-E33 S47-S43
     13 1699  546 (S14W60 L= 21)  1677    298        5 (1636/ 457-1752/ 644) -85.2 W55-W65 S20-S08
     14  574  553 (S46E20 L=300)   253     58        0 ( 554/ 524- 595/ 583) -45.0 E22-E18 S47-S44
     15 1470  476 (S25W46 L=  6)   586     89        0 (1423/ 454-1512/ 502)   0.0 W41-W50 S26-S23
     16 1465  322 (S34W58 L= 19)  2351    283        5 (1325/ 295-1600/ 367)  34.4 W47-W69 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