Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  424 1536 (N13E55 L=254)  4263    779       15 ( 326/1306- 525/1739) -71.0 E63-E50 S02-N27
      2  492 1484 (N11E46 L=263)   814    181        3 ( 446/1408- 541/1548) -77.5 E50-E42 N06-N15
      3  328 1276 (S05E52 L=257)   418    104        1 ( 298/1247- 369/1303) -14.0 E55-E47 S05-S03
      4  134 1180 (S14E72 L=237)  1136    135        0 ( 109/1126- 157/1230)  80.5 E74-E69 S16-S10
      5  463 1071 (S14E36 L=273)   538     93        1 ( 435/1022- 495/1123)  80.5 E37-E33 S16-S10
      6  450  824 (S29E33 L=276)   292     39        0 ( 433/ 794- 472/ 847)  63.4 E33-E30 S30-S26
      7  984  745 (S20W05 L=314)   357     86        3 ( 941/ 719-1035/ 765)  38.7 W02-W07 S22-S18
      8  324  677 (S42E44 L=265)   989    242        2 ( 266/ 632- 389/ 715)  31.6 E51-E37 S46-S37
      9 1709  567 (S11W61 L= 10)  1058    188        4 (1665/ 530-1765/ 594)  29.1 W55-W67 S13-S08
     10 1425  432 (S27W45 L=354)   645    179        2 (1352/ 405-1491/ 466)  10.3 W38-W50 S27-S25
     11 1371  259 (S39W53 L=  2)  1097    238        2 (1340/ 205-1394/ 324) -32.5 W48-W59 S42-S34
  

    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