Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1158 1526 (N34W01 L= 13)  1205    319        3 (1014/1427-1319/1597)  -9.9 E10-W13 N28-N36
      2  753 1542 (N28E28 L=344)   472     96        2 ( 731/1506- 782/1591) -71.6 E29-E25 N25-N31
      3 1603 1474 (N38W37 L= 49)   300    101        1 (1571/1454-1647/1488)   0.0 W33-W40 N36-N39
      4 1375 1410 (N29W18 L= 30)   225     34        1 (1359/1387-1395/1431)  90.0 W16-W19 N27-N30
      5  785 1275 (N11E19 L=354)   516    309        0 ( 689/1193- 892/1383) -25.2 E27-E10 N07-N16
      6  291 1312 (N06E57 L=315)   444     73        1 ( 271/1267- 316/1347)  78.7 E58-E54 N03-N08
      7 1689 1237 (N23W43 L= 55)   955    253        4 (1599/1169-1772/1283) -14.9 W35-W50 N19-N25
      8 1112 1183 (N11W03 L= 15)   619    200        2 (1039/1098-1163/1273)  71.6 W00-W05 N04-N16
  

    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