Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1115 1610 (N35W22 L=130)   661    183        0 (1026/1590-1199/1629) -14.0 W15-W27 N33-N37
      2 1634 1530 (N18W55 L=164)   537     63        0 (1609/1498-1670/1564)  12.5 W52-W61 N15-N18
      3 1724 1531 (N15W65 L=173)   466     36        1 (1699/1512-1743/1550)   0.0 W61-W68 N13-N16
      4 1399 1518 (N23W36 L=145)   267     66        0 (1372/1486-1429/1544)  18.4 W33-W39 N21-N23
      5  429 1310 (N32E35 L= 74)  2010    314        3 ( 308/1219- 548/1381) -31.8 E47-E25 N23-N38
      6  921 1262 (N18E01 L=108)  1316    424        2 ( 722/1209-1082/1373)   2.4 E13-W11 N15-N21
      7  233 1192 (N29E55 L= 54)   482     70        1 ( 211/1164- 258/1219)  18.4 E57-E51 N26-N30
      8 1844  936 (S21W56 L=165)  1780    283        3 (1791/ 814-1888/1089)  69.4 W49-W63 S28-S12
      9 1543  752 (S25W27 L=135)   670    176        0 (1493/ 691-1581/ 834)  41.2 W22-W30 S27-S20
     10  677  615 (S14E30 L= 78)   505    128        2 ( 654/ 546- 693/ 682)  58.0 E33-E27 S17-S09
     11 1194  555 (S29E00 L=108)   778    227        1 (1144/ 460-1244/ 647)  48.0 E05-W04 S34-S23
     12  438  520 (S15E52 L= 57)   724     78        1 ( 413/ 477- 466/ 557)  59.0 E54-E48 S17-S12
  

    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