Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1093 1787 (N43E21 L=172)   874    221        3 (1052/1721-1139/1857) -35.5 E28-E14 N38-N48
      2  533 1603 (N15E48 L=145)   553    168        2 ( 517/1542- 550/1659) -37.9 E52-E43 N11-N18
      3 1093 1395 (N15E06 L=187)   261    205        2 (1056/1369-1127/1450)  80.5 E07-E04 N13-N19
      4 1182 1121 (N02W06 L=199)   381    168        1 (1095/1107-1263/1134)  18.4 W01-W10 N00-N04
      5 1399  924 (S03W24 L=217)  1220    307        4 (1355/ 782-1426/1084) -49.8 W17-W29 S10-N05
      6  224 1026 (S24E55 L=138)   646    222        1 ( 168/ 976- 286/1071)  32.5 E60-E49 S27-S20
      7 1180  811 (S15W14 L=208)  1136    270        3 (1055/ 756-1306/ 883)  49.8 W08-W19 S20-S07
      8  854  535 (S40W04 L=197)  1671    341        4 ( 692/ 507-1033/ 570)  15.6 E09-W16 S44-S35
      9 1446  470 (S25W44 L=238)  1210    251        0 (1323/ 445-1547/ 500)  29.7 W36-W50 S29-S20
     10 1569  414 (S24W57 L=250)   422     46        0 (1545/ 393-1588/ 440)  90.0 W55-W58 S25-S22
  

    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