Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  595 1750 (N27E54 L=143)  3170    531        7 ( 451/1621- 715/1841) -87.1 E65-E42 N15-N35
      2 1009 1772 (N39E26 L=171)   772    135        1 ( 970/1705-1035/1826) -33.7 E32-E20 N34-N42
      3  455 1707 (N21E64 L=133)   604     57        0 ( 430/1681- 483/1743) -45.0 E68-E60 N19-N24
      4 1588 1487 (N37W25 L=222)   420    100        0 (1571/1430-1602/1532) -74.1 W23-W26 N32-N39
      5 1283 1329 (N17W07 L=204)   285    103        3 (1255/1289-1329/1368)  51.3 W05-W09 N13-N19
      6  503 1234 (S07E36 L=160)   232     31        0 ( 488/1216- 524/1259) -45.0 E37-E34 S08-S05
      7 1341 1238 (N13W13 L=209)   232     83        1 (1312/1215-1374/1254)   0.0 W10-W15 N11-N13
      8  134 1182 (S16E73 L=123)  1278    130        3 ( 110/1115- 162/1241)  69.4 E80-E69 S19-S11
      9  267 1138 (S17E54 L=143)   397     78        1 ( 237/1111- 294/1167)  53.1 E56-E51 S18-S14
     10  303 1032 (S22E49 L=148)   391     39        0 ( 284/1004- 321/1064) -76.0 E50-E46 S24-S20
     11  331  991 (S24E46 L=151)   270     53        1 ( 298/ 976- 356/1010)   0.0 E48-E43 S25-S22
     12  221  962 (S28E57 L=139)   351     39        1 ( 200/ 942- 241/ 985)  26.6 E60-E54 S29-S26
     13 1255  930 (S07W15 L=212)   268     51        2 (1229/ 899-1284/ 952)  76.0 W13-W16 S08-S04
     14 1457  734 (S11W33 L=230)   735    208        1 (1362/ 679-1546/ 775)  51.3 W29-W38 S16-S06
     15 1642  633 (S11W51 L=248)   309    102        0 (1622/ 597-1673/ 670) -76.0 W49-W53 S12-S08
     16 1142  601 (S28W19 L=216)   584    256        0 (1017/ 574-1241/ 637)  28.3 W11-W24 S30-S22
     17 1433  450 (S27W46 L=243)  1614    318        3 (1288/ 390-1573/ 492)  15.6 W34-W60 S31-S24
     18 1071  407 (S42W24 L=221)  1087    200        1 ( 979/ 381-1168/ 440)  32.5 W18-W29 S44-S36
     19  522  377 (S63E27 L=170)  1466    221        4 ( 428/ 333- 619/ 420)   1.9 E42-E12 S66-S59
  

    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