Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1082 1648 (N39E10 L=310)  1614    361        3 ( 928/1550-1211/1729) -13.0 E25-W01 N34-N42
      2  510 1581 (N22E49 L=271)  1280    154        3 ( 473/1500- 550/1651) -51.3 E54-E45 N17-N27
      3  252 1482 (N12E71 L=248)  4054    229        6 ( 187/1398- 302/1589) -69.0 E78-E64 N06-N19
      4 1581 1280 (N25W32 L=352)   543    128        0 (1540/1220-1619/1335) -40.6 W28-W35 N21-N27
      5 1159 1284 (N16W03 L=323)   306     64        0 (1130/1256-1178/1328)  78.7 W01-W04 N14-N19
      6  235 1087 (S12E57 L=263)   592     90        2 ( 220/1044- 253/1135) -45.0 E59-E54 S13-S08
      7  194 1055 (S14E61 L=259)   430     79        2 ( 181/1008- 204/1094)  80.5 E63-E59 S17-S11
      8  285 1004 (S16E50 L=270)   911    250        1 ( 210/ 916- 373/1069)  38.7 E58-E42 S22-S10
      9  389 1007 (S14E41 L=279)   279     88        1 ( 361/ 968- 418/1047) -38.7 E43-E38 S15-S11
     10 1619  875 (N01W40 L=  0)   423     67        1 (1581/ 848-1654/ 902)  -9.5 W37-W43 N00-N02
     11 1168  629 (S22W17 L=337)  1653    429        1 ( 990/ 539-1364/ 684)  37.3 W07-W28 S31-S14
  

    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