Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  663 1597 (N36E34 L= 66)   383    121        0 ( 624/1554- 700/1638)  66.8 E35-E32 N32-N39
      2  300 1526 (N27E69 L= 31)  2565    202        8 ( 244/1451- 345/1596) -73.3 E79-E59 N21-N31
      3  382 1536 (N29E58 L= 42)   378     33        2 ( 358/1520- 405/1556)  45.0 E60-E54 N27-N30
      4 1246 1424 (N29W12 L=112)   732    266        1 (1168/1317-1299/1529) -63.4 W06-W16 N21-N35
      5  980 1444 (N28E06 L= 93)  1070    202        6 ( 905/1365-1048/1507) -38.7 E12-E01 N23-N31
      6  451 1360 (N18E43 L= 57)   297     63        1 ( 419/1333- 478/1380)   0.0 E46-E40 N16-N19
      7 1560 1305 (N23W35 L=135)   339     72        2 (1537/1263-1580/1338) -76.0 W33-W37 N20-N25
      8 1830 1137 (N14W60 L=160)   490     84        0 (1811/1094-1852/1172) -90.0 W57-W62 N10-N15
      9  149  908 (S13E70 L= 30)  1228    149        3 ( 124/ 865- 185/ 988) -53.1 E75-E64 S15-S07
     10  185  933 (S11E64 L= 36)   486     87        2 ( 157/ 910- 224/ 968)   6.3 E67-E58 S12-S08
     11  390  924 (S09E42 L= 58)   349    139        0 ( 338/ 880- 448/ 958)  33.7 E46-E37 S12-S06
     12 1197  669 (S19W14 L=114)  2867    724       12 ( 911/ 513-1419/ 897)  42.9 E04-W25 S31-S03
     13  730  800 (S14E17 L= 83)   246     86        0 ( 690/ 779- 770/ 814)   0.0 E19-E14 S15-S13
     14 1726  657 (S17W55 L=155)  1062    258        3 (1656/ 572-1795/ 756)  59.7 W50-W60 S22-S10
     15  313  629 (S29E54 L= 46)  1671    214        2 ( 241/ 543- 390/ 687)  20.1 E74-E44 S36-S24
     16 1725  571 (S22W60 L=160)  1512    181        3 (1651/ 535-1799/ 624)  24.8 W52-W67 S25-S18
     17 1620  526 (S26W50 L=150)   335     43        2 (1603/ 501-1634/ 554) -45.0 W47-W52 S27-S24
     18  604  434 (S41E30 L= 70)  1018    132        3 ( 558/ 392- 670/ 479)  28.3 E35-E22 S44-S36
  

    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