Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  803 1855 (N39E52 L=151)   565     93        0 ( 773/1816- 826/1895) -20.6 E61-E45 N35-N42
      2 1293 1700 (N44E06 L=197)   575    185        0 (1234/1619-1345/1760) -28.3 E12-W01 N39-N46
      3  439 1398 (S00E47 L=156)   829    127        1 ( 412/1333- 471/1447) -32.0 E51-E42 S03-N02
      4 1358 1286 (N18W11 L=214)   752    176        0 (1288/1212-1414/1336) -16.7 W06-W16 N14-N19
      5  936  834 (S19W01 L=203)  1254    467        0 ( 727/ 746-1134/ 878)  43.5 E09-W10 S29-S11
      6  510  710 (S39E25 L=177)  1093    303        1 ( 468/ 581- 582/ 811)  86.4 E28-E21 S46-S30
      7 1042  674 (S25W12 L=214)   577    167        0 ( 969/ 638-1124/ 705)  48.8 W08-W15 S29-S20
      8 1710  651 (S05W56 L=259)   346     44        1 (1687/ 629-1728/ 674)  56.3 W54-W58 S06-S03
      9 1092  552 (S30W20 L=222)   706    224        0 ( 973/ 537-1189/ 570)  31.6 W12-W25 S34-S26
     10  443  427 (S63E35 L=167)  1346    167        4 ( 377/ 390- 516/ 481)  -5.2 E46-E22 S65-S59
     11 1388  430 (S27W45 L=248)   732    191        0 (1320/ 398-1468/ 471)   7.6 W37-W52 S28-S25
     12 1654  414 (S18W69 L=272)   512     78        1 (1620/ 392-1679/ 445)  78.7 W67-W72 S20-S15
  

    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