Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1274 1714 (N47W35 L=134)  1398    286        5 (1239/1618-1325/1783)  38.7 W29-W44 N39-N52
      2 1330 1503 (N31W29 L=128)   374    159        1 (1268/1442-1382/1561) -68.2 W25-W32 N26-N36
      3 1012 1426 (N31W05 L=104)   219     73        0 ( 989/1388-1043/1456) -59.0 W03-W06 N28-N33
      4  742 1377 (N31E15 L= 84)   339    146        0 ( 709/1298- 773/1429) -71.6 E17-E13 N25-N34
      5 1559 1228 (N11W38 L=137)   545    162        4 (1472/1157-1623/1274) -41.6 W32-W41 N05-N14
      6 1265 1172 (N12W17 L=116)   239     80        1 (1226/1143-1295/1210) -51.3 W14-W18 N09-N14
      7 1107 1164 (N13W07 L=106)   264    104        2 (1053/1131-1156/1193)  15.9 W03-W10 N11-N14
      8  238 1046 (N14E57 L= 42)   752    177        2 ( 197/1007- 280/1119) -45.0 E61-E52 N11-N18
      9 1433  548 (S30W22 L=121)  1543    389        1 (1248/ 463-1577/ 674)  17.7 W08-W33 S34-S23
     10  643  474 (S25E34 L= 65)   666    193        1 ( 578/ 440- 707/ 526)  21.0 E40-E27 S27-S21
  

    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