Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1090 1389 (N27W06 L=334)  1234    353        5 ( 917/1330-1242/1436) -21.8 E05-W16 N22-N30
      2 1558 1276 (N17W38 L=  5)   335     80        1 (1545/1233-1573/1321)  68.2 W36-W38 N14-N19
      3 1868 1293 (N14W69 L= 37)   658    120        1 (1838/1264-1898/1313) -18.4 W65-W74 N12-N15
      4 1379 1259 (N17W24 L=352)  1029    188        4 (1316/1203-1447/1309) -40.6 W20-W28 N13-N20
      5  650 1181 (N15E24 L=304)   501    185        2 ( 615/1088- 685/1269) -70.0 E26-E21 N09-N21
      6 1773  877 (S10W53 L= 21)  1635    330        2 (1696/ 784-1845/ 967) -42.9 W45-W61 S16-S03
      7  427  904 (S01E40 L=287)   297     43        2 ( 403/ 883- 451/ 927)  26.6 E42-E38 S02-N00
      8 1391  809 (S11W22 L=350)   356     85        1 (1368/ 765-1419/ 853)  59.0 W20-W24 S13-S08
      9 1486  564 (S28W31 L=359)  2305    409        5 (1301/ 437-1639/ 706)  37.7 W17-W42 S36-S19
     10 1199  624 (S22W10 L=337)   499    165        0 (1123/ 574-1272/ 677)  31.0 W04-W14 S24-S18
  

    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