Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  509 1623 (N29E43 L=320)   559    139        1 ( 491/1566- 534/1678) -63.4 E46-E41 N25-N33
      2  590 1538 (N24E34 L=329)   247     55        0 ( 564/1522- 618/1554)  18.4 E36-E31 N22-N24
      3 1477 1339 (N18W26 L= 29)   568    171        0 (1436/1248-1501/1417) -68.2 W23-W27 N12-N22
      4 1212 1392 (N19W09 L= 12)   290     97        0 (1160/1378-1262/1406)   0.0 W05-W11 N18-N19
      5  897 1310 (N11E10 L=353)   810    269        1 ( 791/1219-1036/1375) -19.4 E17-W00 N06-N14
      6  425 1245 (N04E40 L=323)   343     54        2 ( 407/1214- 448/1278)  90.0 E41-E38 N02-N06
      7 1706 1139 (N09W44 L= 47)   387    134        1 (1662/1078-1744/1186) -45.0 W40-W47 N05-N11
      8  437  724 (S27E38 L=325)  1864    572        3 ( 281/ 585- 613/ 910)  35.4 E55-E24 S36-S14
      9 1479  507 (S32W39 L= 42)   295     73        0 (1439/ 491-1516/ 522)  36.9 W35-W41 S33-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