Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1240 1455 (N23W03 L=221)   313    124        0 (1205/1419-1290/1501) -26.6 W00-W06 N21-N25
      2  546 1425 (N07E38 L=179)   841    178        3 ( 477/1369- 606/1471)  49.4 E41-E34 N02-N09
      3  706 1416 (N09E28 L=190)   452    140        0 ( 644/1366- 764/1447) -11.3 E32-E22 N06-N09
      4 1286 1308 (N15W09 L=227)   589    165        3 (1215/1255-1344/1368) -18.4 W03-W12 N12-N17
      5 1681 1163 (N17W38 L=256)   940    237        2 (1616/1109-1764/1257) -21.0 W31-W46 N14-N21
      6  177 1134 (S16E65 L=153)   900    151        0 ( 153/1051- 202/1191)  66.0 E67-E62 S20-S11
      7  940  812 (S21E00 L=218)   341    169        0 ( 867/ 793-1024/ 834)  36.9 E04-W04 S23-S17
      8  223  733 (S40E61 L=156)   572    143        1 ( 201/ 666- 242/ 796)  56.3 E65-E56 S44-S35
      9  463  661 (S41E33 L=185)   740    118        2 ( 414/ 619- 507/ 700)  37.9 E38-E28 S44-S37
     10 1747  557 (S14W65 L=283)   763     62        2 (1713/ 523-1772/ 585) -63.4 W60-W69 S15-S11
     11 1080  407 (S42W23 L=241)   293     80        1 (1054/ 370-1103/ 442) -23.2 W19-W26 S44-S40
     12 1564  299 (S33W73 L=290)   696     66        0 (1539/ 269-1588/ 329)   0.0 W65-W83 S33-S31
  

    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