Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1129 1649 (N43E04 L=319)   487     81        0 (1102/1609-1153/1687) -45.0 E07-E01 N39-N45
      2  362 1457 (N16E57 L=267)  1473    162        3 ( 331/1383- 404/1527) -78.7 E59-E53 N11-N21
      3  552 1469 (N20E40 L=283)   375     96        2 ( 519/1428- 590/1495) -26.6 E43-E37 N17-N21
      4  865 1393 (N20E16 L=307)   305     97        2 ( 828/1365- 909/1417)  38.7 E18-E13 N17-N21
      5 1693 1330 (N30W45 L=  8)   643    160        1 (1639/1267-1732/1391) -36.9 W40-W48 N26-N33
      6 1039 1271 (N15E03 L=320)   331     94        2 (1008/1238-1076/1320) -38.7 E05-W00 N12-N17
      7 1882 1239 (N27W70 L= 33)   754     75        2 (1869/1192-1899/1275)  50.2 W66-W75 N23-N29
      8 1797 1242 (N25W56 L= 19)   648     44        0 (1775/1209-1819/1269)  63.4 W53-W59 N23-N27
      9 1195  859 (S07W13 L=336)   216    115        1 (1141/ 837-1246/ 884)   0.0 W09-W16 S07-S06
     10 1209  726 (S15W16 L=340)   918    264        1 (1136/ 642-1261/ 780) -83.7 W11-W19 S20-S11
     11 1426  668 (S15W32 L=355)   527    148        0 (1361/ 623-1477/ 712)  49.4 W28-W34 S19-S11
     12  585  638 (S31E24 L=299)  1235    371        3 ( 449/ 585- 760/ 699)  26.6 E37-E11 S37-S24
     13 1314  678 (S17W24 L=348)   258     62        2 (1288/ 642-1337/ 699) -45.0 W22-W25 S18-S15
     14 1133  634 (S22W13 L=337)   222     32        1 (1112/ 614-1153/ 660) -45.0 W12-W15 S23-S20
     15 1543  444 (S28W52 L= 15)  3770    478        5 (1346/ 375-1732/ 553)  32.2 W37-W65 S35-S17
     16  989  380 (S43W10 L=334)  1450    302        2 ( 872/ 349-1129/ 408)  20.2 W00-W20 S46-S38
  

    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