Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1642 1449 (N20W52 L=186)   570    148        1 (1592/1425-1692/1469) -20.6 W46-W55 N17-N21
      2  860 1217 (N20E07 L=127)   284    132        1 ( 819/1153- 897/1285)  45.0 E10-E03 N16-N24
      3 1164 1191 (N14W11 L=146)   481    135        0 (1139/1116-1186/1257)  83.7 W09-W12 N09-N18
      4 1605 1026 (S04W37 L=172)   548    135        1 (1547/ 996-1661/1050) -20.6 W33-W41 S05-S01
      5 1893  890 (S19W63 L=198)  1863    242        3 (1861/ 808-1935/ 958)  32.5 W59-W70 S24-S14
      6 1666  826 (S17W39 L=174)   980    224        1 (1612/ 718-1693/ 930)  59.7 W34-W42 S23-S10
      7 1286  676 (S19W11 L=146)   378    108        2 (1262/ 624-1303/ 732)  66.8 W09-W12 S22-S15
      8  297  689 (S06E58 L= 77)   363     61        1 ( 280/ 667- 324/ 709)   0.0 E59-E54 S06-S04
      9 1577  634 (S28W32 L=166)   342     76        1 (1551/ 606-1602/ 657)  26.6 W29-W33 S29-S25
     10  342  608 (S11E56 L= 79)   463     54        2 ( 322/ 582- 364/ 636)  71.6 E57-E53 S12-S09
     11  780  506 (S22E25 L=110)  1073    256        5 ( 680/ 451- 879/ 564)  14.0 E33-E17 S25-S19
     12  565  320 (S34E50 L= 85)   811    138        0 ( 534/ 292- 589/ 357)  32.0 E53-E45 S36-S30
  

    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