Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  830 1648 (N39E25 L=180)   432    101        0 ( 789/1623- 875/1669)  -7.1 E29-E20 N37-N40
      2 1588 1481 (N37W39 L=245)  2022    348        3 (1543/1340-1620/1624)  79.2 W35-W45 N27-N48
      3 1074 1511 (N33E03 L=202)   669    227        3 (1025/1415-1146/1593) -47.7 E07-W03 N26-N38
      4  589 1480 (N24E37 L=168)   425    151        0 ( 557/1435- 645/1504)  38.7 E39-E33 N21-N26
      5  915 1401 (N23E12 L=193)   294     85        1 ( 879/1377- 956/1423)  -9.5 E14-E08 N21-N24
      6 1250 1211 (N15W12 L=217)   278    119        0 (1221/1169-1301/1257) -45.0 W10-W15 N12-N17
      7  616 1123 (N02E27 L=179)   382    104        2 ( 596/1077- 635/1178) -68.2 E28-E25 N00-N05
      8  364  928 (S13E44 L=161)   639    103        1 ( 329/ 890- 395/ 958) -36.9 E47-E41 S14-S10
      9 1524  886 (S02W34 L=239)   381     92        1 (1479/ 851-1557/ 918) -18.4 W30-W36 S03-N00
     10 1003  841 (S10W01 L=206)   211     41        0 ( 977/ 819-1023/ 864)  63.4 W00-W01 S12-S08
     11 1319  742 (S13W22 L=227)   850    232        1 (1243/ 685-1358/ 818)  60.9 W17-W24 S17-S07
     12 1235  771 (S12W16 L=221)   217     42        0 (1207/ 753-1259/ 788)  45.0 W14-W17 S13-S10
     13  532  695 (S26E30 L=175)   248     46        0 ( 514/ 678- 557/ 721) -33.7 E31-E28 S26-S23
     14 1192  586 (S24W17 L=222)   490    104        2 (1138/ 552-1229/ 614)  -9.5 W13-W19 S26-S22
     15  458  525 (S39E40 L=165)   327     27        1 ( 440/ 505- 475/ 549)  71.6 E41-E38 S40-S37
     16  731  535 (S34E15 L=190)   256     66        0 ( 693/ 522- 765/ 547)  11.3 E17-E12 S35-S33
     17  375  494 (S43E53 L=153)  1550    145        3 ( 299/ 470- 432/ 536)   7.4 E68-E44 S45-S39
     18 1048  234 (S55W19 L=224)   362     36        1 (1026/ 214-1069/ 252)  45.0 W15-W20 S57-S52
  

    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