Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  533 1676 (N26E53 L=175)   688    182        0 ( 488/1625- 562/1706) -21.8 E58-E48 N22-N27
      2  846 1433 (N17E21 L=208)   233     89        1 ( 833/1381- 861/1471) -59.0 E22-E18 N13-N19
      3 1190 1354 (N20W02 L=231)   436    169        0 (1178/1267-1204/1430) -74.7 W00-W04 N14-N25
      4 1475 1314 (N25W22 L=250)   512     74        1 (1438/1276-1508/1349) -18.4 W18-W24 N22-N26
      5 1101 1165 (N07W01 L=230)   561    215        2 (1046/1060-1140/1215) -45.0 E02-W05 N01-N09
      6 1615 1015 (N11W37 L=265)   627    154        3 (1562/ 969-1666/1076) -11.3 W31-W41 N08-N13
      7  446  846 (S26E34 L=195)   789    272        2 ( 391/ 797- 491/ 926)  66.0 E38-E30 S29-S20
      8 1568  584 (S15W47 L=276)   699    101        3 (1521/ 553-1614/ 610)  33.7 W43-W49 S16-S12
      9  936  565 (S32W06 L=235)   346     78        1 ( 911/ 542- 977/ 601) -21.8 W04-W09 S34-S29
     10 1430  402 (S29W46 L=275)   704    106        3 (1395/ 368-1470/ 434)  45.0 W43-W49 S31-S26
     11 1161  299 (S44W35 L=263)   410     80        0 (1139/ 272-1179/ 326) -18.4 W31-W37 S45-S42
     12 1212  253 (S46W42 L=271)   597    107        0 (1174/ 222-1248/ 295) -81.9 W40-W44 S49-S41
  

    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