Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  561 1677 (N20E50 L=223)   854    124        4 ( 504/1649- 616/1715)  29.7 E53-E45 N17-N22
      2  866 1452 (N13E20 L=253)   664    310        0 ( 818/1359- 913/1601) -42.9 E28-E14 N08-N21
      3  200 1192 (S15E61 L=212)  4051    848       20 ( 121/ 999- 290/1349)  -4.9 E69-E52 S27-S05
      4  696  990 (S17E18 L=255)   686    218        0 ( 599/ 952- 810/1017)  12.1 E24-E10 S19-S15
      5  169  942 (S31E59 L=214)   714    127        2 ( 134/ 922- 205/ 960)   0.0 E63-E55 S31-S29
      6  387  939 (S27E39 L=235)   486    142        2 ( 333/ 921- 442/ 957)  12.5 E43-E34 S28-S24
      7 1371  872 (S06W23 L=297)   605    230        0 (1295/ 790-1424/ 925) -14.0 W17-W29 S08-S04
      8  919  804 (S22E00 L=273)  1370    555        2 ( 642/ 706-1156/ 885)  15.3 E17-W16 S30-S16
      9 1367  593 (S21W33 L=307)   332     90        0 (1332/ 554-1394/ 626)  -9.5 W30-W36 S22-S19
     10 1098  465 (S36W23 L=297)   508    105        2 (1046/ 422-1140/ 510) -18.4 W17-W27 S38-S35
  

    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