Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  971 1815 (N52E26 L=199)   426     37        1 ( 952/1793- 995/1835) -18.4 E29-E23 N49-N54
      2 1038 1725 (N46E15 L=211)  1182    147        4 ( 980/1684-1100/1775)  11.3 E18-E08 N42-N50
      3 1047 1605 (N36E10 L=216)   568     87        1 (1007/1577-1091/1639)  31.0 E11-E06 N33-N39
      4  371 1533 (N20E59 L=166)  1269    246        1 ( 299/1446- 456/1617) -14.4 E81-E46 N15-N24
      5  275 1377 (N09E61 L=165)  1279    147        3 ( 230/1302- 316/1431)  77.5 E64-E56 N03-N12
      6 1331 1178 (N13W16 L=242)   626    133        2 (1268/1132-1389/1218)  45.0 W12-W19 N09-N16
      7 1277 1110 (N09W14 L=239)   884    281        2 (1209/1047-1381/1213) -22.6 W08-W20 N05-N13
      8  769 1173 (N04E18 L=208)   235     75        0 ( 736/1154- 805/1197) -21.8 E20-E15 N03-N05
      9 1058 1114 (N05W01 L=226)   435    122        0 (1038/1057-1079/1181) -74.1 E01-W02 N02-N09
     10  397 1003 (S12E41 L=185)  1040    250        0 ( 378/ 877- 420/1123) -78.7 E42-E37 S19-S04
     11 1349 1017 (N04W20 L=245)   219     64        0 (1324/ 993-1368/1043)  63.4 W18-W20 N02-N06
     12 1554 1011 (N07W33 L=259)   420    130        0 (1504/ 969-1596/1038)  -8.1 W29-W36 N05-N08
     13 1290  749 (S13W21 L=246)   950    342        2 (1221/ 635-1374/ 914)  84.0 W18-W22 S20-S01
     14  306  856 (S22E48 L=178)  1060    188        0 ( 221/ 816- 391/ 877)  15.5 E58-E40 S25-S19
     15 1672  604 (S15W54 L=279)   333     56        0 (1644/ 576-1693/ 620)  -9.5 W50-W56 S16-S13
     16 1503  419 (S30W49 L=275)   505    112        0 (1460/ 392-1552/ 442)  39.8 W45-W52 S32-S27
     17  487  421 (S50E37 L=188)   616     97        0 ( 442/ 406- 536/ 439)  12.1 E44-E30 S52-S48
     18  656  336 (S54E19 L=207)  1367    243        2 ( 571/ 271- 734/ 421) -38.2 E27-E10 S59-S47
     19  964  237 (S56W16 L=242)  1047    121        2 ( 916/ 188-1021/ 291) -26.6 W09-W25 S59-S51
  

    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