Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  656 1711 (N29E46 L=139)  1062    148        3 ( 599/1662- 710/1742)  29.7 E50-E42 N25-N31
      2  375 1580 (N14E62 L=123)  2856    332        5 ( 321/1500- 421/1660)  90.0 E69-E56 N08-N19
      3 1085 1603 (N33E13 L=173)   276     78        0 (1049/1572-1129/1618)  -8.1 E15-E08 N31-N33
      4 1073  839 (S13W07 L=193)   506    156        4 (1015/ 786-1121/ 896)  74.1 W04-W10 S16-S09
      5 1634  749 (S03W45 L=231)   465     63        1 (1609/ 713-1654/ 779)  90.0 W43-W47 S05-S01
      6 1291  720 (S14W24 L=209)   447    126        0 (1228/ 652-1346/ 776)  73.3 W21-W25 S19-S09
      7  750  554 (S39E05 L=180)  1156    330        3 ( 606/ 526- 919/ 581)  17.7 E16-W06 S42-S33
      8 1469  522 (S20W43 L=229)   362    135        0 (1426/ 504-1512/ 549)  51.3 W39-W45 S22-S17
      9  318  523 (S54E49 L=137)   590    167        0 ( 278/ 507- 364/ 543)  14.0 E54-E42 S55-S51
     10 1237  415 (S33W33 L=219)   486    147        1 (1186/ 388-1316/ 444)  26.6 W28-W37 S35-S29
     11 1572  390 (S25W59 L=244)  1869    446        2 (1465/ 334-1652/ 424)   7.8 W49-W71 S27-S22
  

    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