Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  751 1642 (N52E11 L=249)   528     77        3 ( 719/1607- 791/1665) -20.6 E14-E06 N49-N54
      2 1824 1394 (N12W72 L=332)  1096    151        3 (1780/1350-1854/1441)  20.0 W66-W80 N09-N15
      3 1761 1300 (N08W58 L=318)   811    151        4 (1706/1258-1811/1351) -45.0 W54-W63 N05-N12
      4  201 1088 (N20E63 L=197)  1390    215        5 ( 135/1048- 289/1129)   0.0 E74-E52 N17-N22
      5 1167  887 (S04W06 L=266)   488    220        3 (1076/ 860-1270/ 901) -20.0 W01-W12 S05-S01
      6 1546  795 (S17W30 L=290)   240    110        2 (1535/ 745-1561/ 852) -90.0 W29-W31 S20-S13
      7 1280  501 (S31W08 L=268)  3175    648       10 (1007/ 360-1475/ 688)  20.7 E14-W23 S37-S21
      8  268  618 (S10E65 L=194)   524     59        2 ( 252/ 595- 293/ 641) -26.6 E68-E62 S11-S08
      9  917  571 (S20E14 L=245)   219     84        0 ( 887/ 538- 945/ 599) -68.2 E15-E13 S22-S17
     10  518  548 (S16E43 L=217)   365     96        2 ( 494/ 506- 552/ 576)  29.7 E46-E39 S18-S14
     11  737  411 (S28E32 L=228)   608    180        2 ( 718/ 344- 763/ 502)  67.4 E35-E27 S33-S21
     12  523  349 (S31E55 L=205)   746    123        1 ( 501/ 290- 544/ 405)  35.5 E63-E48 S36-S26
     13  803  270 (S41E35 L=225)   409     58        1 ( 773/ 249- 834/ 292)   8.1 E38-E31 S42-S39
  

    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